arXiv: 1504.07123v2 [quant-ph] 8 Oct 2015 


Nonclassical properties and quantum resources of hierarchical 
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We motivate and introduce a class of “hierarchical” quantum superposition states 
of N coupled quantum oscillators. Unlike other well-known multimode photonic 
Schrodinger cat states such as entangled coherent states, the hierarchical superpo¬ 
sition states are characterized as two-branch superpositions of tensor products of 
single-mode Schrodinger cat states. In addition to analyzing the photon statistics 
and quasiprobability distributions of prominent examples of these nonclassical states, 
we consider their usefulness for high-precision quantum metrology of nonlinear opti¬ 
cal Hamiltonians and quantify their mode entanglement. We propose two methods 
for generating hierarchical superpositions in = 2 coupled microwave cavities which 
exploit currently existing quantum optical technology for generating entanglement 
between spatially separated electromagnetic field modes. 


I. INTRODUCTION 


Maximally entangled quantum states occupy a distinguished position in the theory of 
quantum information. One has only to consider the central role of Greenberger-Horne- 
Zeilinger (GHZ) states of Q in many quantum algorithms and quantum teleportation 

protocols to be convinced of their importance. From a practical perspective, such states 
comprise the most valuable quantum resource, both in terms of entanglement and usefulness 
for quantum metrology js, 4], Unfortunately, in the case of N two-level quantum systems 
such as spin-1/2 chains, the maximally entangled states are sensitive to local errors (e.g., 
phase flips) and can quickly lose all nonclassical resources. However, because of the countably 
infinite dimension of the Hilbert space of a chain of quantum oscillators (isomorphic to 
(C(C))®^), one may hope to engineer maximally entangled states of a subspace isomorphic 
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to that are robust under quantum evolutions corresponding to relevant sources of 

decoherence. 

The entangled coherent states (see Ref.js, ^ and references therein) are paradigmatic 
examples of entangled states of £^(C) 0 £^(C). However, only a strict subset of these are 
maximally entangled. Explicit conditions for maximal entanglement of linear combinations 


of products of coherent states have been found and multimode entangled coherent states 
have been studied in the context of encoding continuous variable quantum information |S- 
Q. To extend the notion of GHZ states to separable Hilbert space (i.e., Hilbert space with 
countable orthonormal basis) in a general way, we first introduce the following two-branch, 
N mode states: 

+ ( 1 ) 


1 


^2 + 2Re(z^) 


where |0) is a single-mode pure quantum state of Hilbert space Ti, and [/ is a partial 
isometry with 10) in its domain. When z = 0, the branches 10)®-^ and 17®'^|0)®'^ are 
orthogonal in and the resulting state is a strict analog of a GHZ state. For example, 
superpositions of oscillator Fock states proportional to , with ip e [0, 27r], fall 

into this class and are interesting for their quantum optical properties. Such superpositions 
represent the most obvious generalization of GHZ states to the Hilbert space of N oscillators, 
(£^(C))®'^. The case of z 0, although deviating from the strict notion of GHZ states 
due to non-orthogonality of the branches, contains many important macroscopic Wmode 
superpositions. The entangled coherent states (having |0) = |q;) and U = e*^“^“T>(/3), where 
0 G [0, 27r] and T>(/3) is the oscillator displacement operator for a, /3 & C) serve as well- 
studied examples. However, the focus of this paper will be on revealing certain new states of 
the z = £) set. In particular, consider taking |0) = |ei) := ^ where V 

a single-mode partial isometry containing |0') in its domain and w = (0'|F|0') G M. Then 
the state |e 2 } := (^ j 2 - 2 w } orthogonal to |ei). Taking U = e*®/'^(|ei)(e 2 | + |e 2 )(ei|) 

produces two-branch superpositions of the following form, which we refer to as hierarchical 
cat states (HCS): 


|HCS^) 


_L ( /(i+yMT" I 

V2\\ V2T^ ) \ ) ) 


( 2 ) 


The origin of the name “hierarchical cat state” is self-evident: |HGS^) is an equal weight 
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superposition of two orthogonal branches in (i.e., an A^-mode “cat” state), while each 
branch is a tensor product of “kitten” superpositions in the single-mode Hilbert space T-L. 
This construction allows to consider maximally entangled states of a subspace of 

(£^(C))®'^ that retain single-mode quantum coherence in the states |0') ± V\(p') even after 
intermode decoherence processes reduce the superposition state |HCS^) to an equal proba¬ 
bility statistical mixture of tensor products. 

When l-L = ^^(C), |HCS^) is in a two (complex) dimensional subspace of (£^(C))®'^. 
By appropriate choices of V and |0'), the branches of |HCS^) can take the form of tensor 


products of single-mode pho 


states \'ip±) oc |q;) ± 


-a) |n 


jonic Schrodinger cat states such as the even and odd coherent 


mi or superpositions of squeezed states. For these photonic HCS 
states, it is clear that if the single-mode coherence time is sufficiently long (e.g., greater 
than the intermode coherence time), a statistical mixture of Wmode Schrodinger cat states 
remains even after intermode coherence is lost by some decohering process. 

The photonic HCS state obtained by taking V = to be the oscillator tt phase 

shift, 10') = la) a coherent state of the quantum oscillator with amplitude |a|, and d = 0 
or 6* = TT in (Ej) was introduced in Ref.|l2|. We will label these states by |HCS^(a)) (the 
-|- symbolizing 0 = 0 and the — symbolizing 6 = tv) and they will serve as the canonical 
examples of photonic HCS in this work. Each branch of |HCS^(a)) is an Wfold tensor 
product of either the even coherent state 10+) or odd coherent state |0-). If A/” is an odd 
natural number, then the even (odd) branch is an eigenvector of the photon parity operator 
-with eigenvalue 1 (—1) and hence, “^“0HCS^(a)) = |HCS^(a)). If N is 

an even natural number, |HCS^(a)) is invariant under such a local tt rotation. Independent 
of N, a photon parity measurement results in a projection of |HCS^(a)) onto either the even 
or odd branch. The entangled coherent states |ECS^(a)) oc |a)®-^ ± |—a)®-^ are invariant 
under the bosonic algebra freely generated by OjUj, for i, j G |1,.. • A^}, so that they are 
considered to be Barut-Girardello coherent states of 5p(A^, C) ISj. The state |HCS^(a)) is 


invariant under the algebra freely generated by the identity operator and the two-photon 
annihilation operators a|, j = This algebraic property allows |HCS^(a)) to be 

considered as superpositions of sp{N,C) Barut-Girardello coherent states. However, while 
such properties as quasiprobability densities, photon statistics, and dissipative evolutions of 
the entangled coherent states have been thoroughly documented 
of the properties of photonic HCS states is lacking. 


1f)| , detailed analysis 
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When |HCS^(«)) is shared between two spatially separated parties, the state |HCS^(a;)) 
serves as an entanglement resource for teleportation of an arbitrary superposition of coherent 
states of the form ci|q;) + C 2 I—a) in the same way that the GHZ state 

serves as an entanglement resource for teleportation of an arbitrary qubit pure state. This 
follows from the fact that |HCS^(q;)) are maximally entangled states in the 4-dimensional 


sub-Hilbert space spanned by {|(—lya) ® |(—The states |HCS 


a 


) are 


also useful probes for high-precision phase estimation of Hamiltonians of 5((2,C) (l2| (see 
also Section IHT At . These intriguing attributes motivate a more thorough description and 
analysis of HCS states. 

The remainder of this paper is structured as follows: in Section El we indicate some 
nonclassical properties of |HCS^(q;)), citing the nonclassical properties of the entangled co¬ 
herent states for comparison; Section HI is devoted to exploration of the quantum resources, 
in particular the metrological usefulness and entanglement entropy, of |HCS^(a;)); in Section 
IV, we propose two methods for generating |HCS^(q;)) using techniques which are acces¬ 
sible by current quantum optical technology; in Section V, we demonstrate how the idea 
of hierarchically encoding continuous variable quantum information can be deepened with 
many levels of hierarchy. We do not attempt an exhaustive analysis of photonic hierarhical 
cat states, but rather try to show the most salient properties of these states by considering 
basic examples. 


II. NONCLASSICAL PROPERTIES OF HCS 


Here, we take note of the basic photon statistics of |HCS^(a;)) for arbitrary N, derive some 
of its quasiprobability distributions for V = 2, and show a duality between Pauli matrices 
and photon operations in the subspace of £^(C) spanned by the even and odd coherent 
states. Throughout this section, we compare the nonclassical properties of |HCS^(q;)) to 
those of the entangled coherent states, which are more familiar. We show the inner products 


^ In fact, the states |HCS^(a)) exhibit the same amount (1 ebit) of entanglement entropy as IECS 2 (a)) 


16, 


17|. A “Bell basis” of maximally entangled states for this 4-dimensional sub-Hilbert space is given 


by: |HCS+(a)), |HCS 2 -(a)), |ECS 2 -(a)), ® I)|ECS 2 -(a)) [l 
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of |ECS^(«)) and |HCS^(a;)) immediately: 
{ECS+(a)|HCS+(a)) = ,, ^ 


{ECS]!;(a)|HCS;;(a)) = 


1 


Vl + e-' 


2Ncfl 


(ECS^(a)|HCS±(«)) = 


(1 + + (1 

_ lg-2a2^Ar/2 

,N even 

(3) 

(i + 


, N odd 




, N even 

(4) 


_ lg-2a2^Ar/2 

,A^ odd 


f 0 

N even 




,A^ odd 


( 5 ) 

\ \/ 1 —e“2Aro;2 





where we have taken a e M. It is intriguing to take note of the a —>■ oo asymptotics of these 
inner products (the N ^ oo asymptotic is always zero). Eor any odd N and for reasonably 
large a, (ECS)()(a)|HCS)i/'(«)) ~ 0, while for any even N and any a, (ECS)^(a)|HCS^(a)) = 
0 identically. The total expected photon number in all of these states is asymptotically Na"^, 
i.e., |HCS^(a)) and |ECS^(a)) differ mainly in photon statistics and quantum correlations, 
not in intensity. In particular, if Pg (Po = I — Pe) is the projection onto the even (odd) 
photon number subspace^ of (£^(C))®'^, it is clear that Po|ECSj^(a)) = Pe|ECSj^(a)) = 0 
while P<,|HCSj(,(a)) = Pe|HCSj^(a)) = 1 /^ 2 . 

Because of its symmetry, |HCS^(a)) has a simple expression as a superposition of tensor 
products of coherent states. To do this, first consider the direct product Z 2 x • • • x Z 2 {N 
times) with group operation addition modulo 2. This group is isomorphic to the abelian 
group It with elements given by , kj G {0,1}, and group operation being 

operator multiplication. Let Hi (II 2 ) be the subgroup of elements corresponding to k such 
that the number of nonzero entries of k is even (odd). Then, one can write (again for a G M): 

^Na^/2 ( , . \ 

|HCS^((a)) = I ^cosh“^/^ + (—sinh“^/^ ^ 1 (^) 

From the above expression, the expansion in the Fock state basis can be made explicit by 
making use of the fact that \a) = 


Explicitly, Pg = X]n ||n||2=o mod 2 whcrc Pfi is the rank one projector onto the ray |ni) (8) • • • (8) |nAr). 
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P{n,m) 



FIG. 1: Photon number distribution P{n,m) for |HCS]!^^ 2 ('^ ~ ^))> fn ^ {1,. ■ ■ ,20}. 


The photon number distribution of |HCS]!}(q;)) for any a G C is given by: 


Pin) = |(n|HCS+(a)}p 

„-V|a|2 


2N+1 


N 


E 




JJ(1 + (_!)...+,) 


a 


rik 


^riki 


( 7 ) 


where |n) := |ni) 0 • • • 0 |niv) for n G N-^. It is clear from the above expression that if n 
does not have all entries even or odd, then the photon number distribution vanishes. For 
N = 2, this results in a checkerboard pattern of zero and nonzero probabilities on the lattice 
Z>o X Z>o (FigH]). This feature stands in contrast to the distribution |(n|ECS^(Q;))p, which 
is identically zero if and only if is odd. Our present focus on the photon statistics 

of |HCS^(q;)) is merely due to the fact that they are the hierarchical cat states of most 
immediate practical use for continuous variable quantum information processing. Indeed, 
more complex photon statistics are furnished by hierarchical cat states formed from, e.g., 
the Z/4Z coherent states, which generalize the even/odd coherent states by forming a 
subspace of £^(C) having orthonormal basis 
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FIG. 2: Mandel Qm parameter for |HCS^(q!)) for |q!| € [0,3]. 
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( 8 ) 


where {n\ej) ^ 0 if and only if n = j mod 4. 

The variance of a single mode quadrature + a]e*®) in |HCS^(q;)) is 


1 

2 


+ |Q;p(coth 2 |q;P + cos(2ArgQ; — 26*)) 


( 9 ) 


for all 6*, i.e., since this variance is greater than 1/2, |HCS^(a)) is not squeezed in any 
quadrature. |HCS^(«)) do not exhibit second order squeezing 2Qj (also referred to as am¬ 
plitude squared squeezing) due to their being eigenvectors of for all j = 1,... ,N. The 
single mode Mandel Qm parameter of |HCS^(q;)) is always negative, but exhibits a dip 
for I a I ~ 3/2 as shown in Figl2l). We now turn to the quasiprobability distributions and 


functional representations of [HCS^ («)) with the focus being on the functional form of the 













hierarchical coherences. In addition, we utilize the explicit expressions of the quasiprobabil¬ 
ity distributions to infer nonclassical features of the states. Each quasiprobability distribu¬ 
tion for a quantum state p of (t'^(C))®^ is obtained by Fourier transformation on of a 
quantum characteristic function corresponding to a particular ordering of bosonic operators 
2 l| . It is well known that negative values of the singular quasidistribution (i.e., Sudarshan- 
Glauber function) for a given quantum state p indicate nonclassical photon statistics, i.e., 
indicate that the photon number distribution is not Poissonian j^]. In a similar spirit, the 
existence of negative values of the Wigner function for a given state indicates non-Gaussian 
quadrature correlations. The explicit form for the Wigner function of |HGS)!)^(q;)) was shown 
in Ref. jl^ . but the analytic expression is only useful for technical purposes. The important 
point is that the Wigner function of a photonic state has an interpretation as a continuous 
set of interference experiments. This is clear from the definition of the single mode Wigner 
function as W{'^) = ( 2 / 7 r^)(D( 7 )e*™^“D(— 7 )); it shows how a quantum state changes when 
its coherent state components are displaced by —7 in phase space, then are reflected in 
phase space, and are displaced again opposite to — 7 . Goherence between coherent state 
components appears as fringes in the Wigner function. Accordingly, the states |HCS^(q;)) 
exhibit two sources of interference: 1 ) from the coherence in each branch and 2 ) from 

the coherence between the branches. These two sources of coherence are not immediately 
visible from the expression of IHGS 2 («)) as an unequal superposition of tensor products of 
coherent states. 

The existence of zeroes of the Q-function of a given quantum state indicates that the sin¬ 
gular quasidistribution for this state takes negative values and, hence, exhibits nonclassical 
features. The Q-function of |HGS)^(a;)), which is a true probability distribution on C^, is 
given by (5(|HGS)^(«))) = |HGS)^(q;))) p. For N = 2 and a G C we have: 


TT 


^|HCS+(a))(/^l>/^2) — W (Q+(/5 i)Q+(/52) + Q-(/3i)Q-(/52)) 


=-(|/3i|W|/32|h 


+ 


27r sinh(2|Q;P) 

— sin(2Im(/5iQ;)) sin(2Im(/d2Q;))) 


(sinh( 2 Re(/diQ;)) sinh( 2 Re(/ 52 Q;)) 


( 10 ) 


where ( 5 +(-)(^) is the Q-function of the even (odd) coherent state . The Q-function 
of |HGS^(a;)) vanishes if and only if each of the terms vanishes. Let us take (3i, (^2 G 
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C such that Re(/3i«) = Re(/32Q;) = 0. In addition, we require that /3i, /32 satisfy: 1) 
Im(/5iQ;) = where k E Z, and 2) Im(/d 2 Q;) = mn, where m E Z. Under these 

constraints, Q+(/di) = Q-{j32) = 0. For these valnes, sinh(2Re(/3iQ;)) sinh(2Re(/52Q;)) — 
sin(2Im(/5iQ;)) sin(2Im(/d2Q;)) = 0, and hence the Q-fnnction is zero at these points of C x C. 

In addition to the quasiprobability distributions, the qnadratnre distribution, calculated 
as the square modulus of the Schrodinger wavefunction, is an especially useful true proba¬ 
bility distribution for systems of oscillators. However, the quadrature distribution is quite 


specific; all that is needed is a representation of the pure state IHCS 2 («)) in a functional 
Hilbert space. We choose the Bargmann representation [231] becanse of the fact that snch 
relevant qnantities as the Schrodinger wavefnnction and the singnlar qnasidistribntion of 
a pure state can be derived from the Bargmann representation by the use of variants of 
the Segal-Bargmann transformation. As an analytic function /|Hcs*(a)} on C x C, the state 
IHCS 2 («)) is represented by 

f , (, ( cosh(a:(^ + w)) qp cosh(a:(^ - w)) \ 

/|HCS±(a))U>«'l - V^e I ;^_g-4|a|2 I' 


That /|HcsJ(a)>(^>'^) f^’kes the form of a sum of unequally weighted functions is a conse¬ 
quence of the fact that IHCS 2 (a)) is an nneqnally weighted superposition of tensor prodncts 
of coherent states. 

The subspace /C = of P(C) spanned by the even and odd coherent states \'ip±) has 
the property that certain photonic operators carry ont equivalent operations as the Pauli 
matrices in this subspace. This allows for quantum operations of a two-level system to be 
interpreted as photonic operations compressed to this subspace. For example, keeping in 
mind the action of \ + |'0_)('^+| in JC we note that a|'0+) = ay/tanh |Q;p|'0_) = 

a-y/tanh \a\‘^ax\'ip+)- Considering the Pauli matrices as observables of a spin-1/2 particle, 
we find the following expressions in terms of self-adjoint photonic observables: 




(^Z 


\a\ 


fl!lV!!ThhElp^j,(i+Arg«)j>^ 

\a\ 


( 12 ) 
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where P/c is the projection to fC. These expressions for Pauli matrices show a duality between 
quadratures and “magnetization” in the subspace /C. For example, if |q;| is sufficiently 
large, we have P/c = V^jafax, i.e., the interaction picture dynamics of a quantum 

oscillator (restricted to /C) driven with unit amplitude is equivalent to a spin-1/2 particle 
with magnetic moment v^lctl in a unit magnetic held along the x-axis. The Pauli matrices in 
(TT^ do not have unique expressions in terms of products of photonic operators and P/c, due 
to the fact that one can rewrite P/c as (l/2|a;p)P,ic(a^e“^*^'^®" + Some of these 

alternative expressions can be instructive; for example, one can rewrite = Pjc cos{T:a^a)Pic- 
A similar duality can be derived for 5u(4) observables in terms of projectors in the subspace 
spanned by the list (jH|) and photonic operations. 


III. QUANTUM RESOURCES OF HCS 


A. Metrological usefulness 


We begin this section by recalling the main problem of quantum metrology and how 
certain quantum states can be utilized for estimation of dynamical parameters at higher 
precision than any classical states. Given a smooth manifold M, let quantum states be 
encoded by a differentiable map specihed by A pa for A G M. The goal is to estimate 
the parameter A with greatest possible precision by making use of an optimal quantum 
measurement and optimal classical post-processing of the measurement results. 

In the present section, we are concerned with the special case of estimation of a displace¬ 
ment parameter. In this case, the parameter manifold is a line with real coordinate x G M 
and the state := lies on a path parametrized by x and generated by the 

self-adjoint, x-independent operator H. If {M{dx)} is a positive operator-valued measure 
(we will refer to {M{dx)} as a “quantum measurement” or, simply, “measurement” from 
now on) which is an unbiased estimator of x, i.e., x = J^x'tv(pxM{dx')), then the quantum 


Cramer-Rao theorem 


24, 


2f)l | states 


{{Sxuf) > 


(13) 


tr(poT2) 

where {(SxmY) ■= fj^(x' — x)^tr(pxM(dx')) is a general expression for the variance of the 
quantum measurement and where L = is the symmetric logarithmic derivative operator 
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defined by the equation ^ = —i[H,p] = \{Lp + pL). The quantity tr(po-h^) is called the 


26|. 


dx ^ h'j 2 ' 

quantum Fisher information of po and is constant on the unitary path generated by H 
If Po = |'0o)(t^o| is pure, then tr(poT^) = 4(tr(po-f^^) — tr(po-f^)^)- Hence, if an experimenter 
has unconstrained access to measurements saturating inequality flTHl) , then a quantum state 
with larger quantum Fisher information with respect to H can be considered as a more use¬ 
ful resource for estimating the displacement parameter x. In this section, we focus on using 
certain multimode pure states l-^o) £ as probes for displacement metrology for paths 

generated by 1-local Hamiltonians , i.e., H having the form H = X]j=i ® Iat-i- Specif¬ 
ically, the metrological problem at hand consists of: 1) preparation of N oscillator modes 
in the probe state |ECS^(a;)) or |HCS^(«)), 2) application of a global unitary operator 
with Hj = H &IV oscillator Hamiltonian and x G M, and 3) estimation of x by an 
optimal separable measurement on the N modes. It is important to note that determination 
of the optimal separable measurement corresponding to the probe state and Hamiltonian 
H requires methods of quantum estimation theory; in particular, the optimal measurement 
does not necessarily correspond with traditional methods of oscillator signal detection such 
has homodyne detection. In this section, we assume that the optimal measurement can be 
performed for any H and we determine the set of H for which |ECS^(«)) and |HCS]!^(a;)) 
allow for a greater precision in the determination of x than the tensor product branch states 
|±q;)®'^ and , respectively. 

As an example of displacement estimation in a finite dimensional Hilbert space, one can 
consider the problem of estimation of a phase parameter 9 imprinted on a quantum state 
P 0 = poe"^^^. Let us take H = ® Iat-i and take po to correspond to the 

GHZ state lO)®-^ -|- |l)®'^/v^. The quantum Eisher information of po with respect to H 
is 4A^; in fact, this is the maximum possible value of the quantum Fisher information 


in with respect to 1-local Hamiltonians of unit operator norm [3|. In contrast, 

any product state has maximal quantum Fisher information of order N over the set of 
such Hamiltonians. This fact suggests an ordering of superposition states based on their 
maximal usefulness for quantum metrology as compared to the maximal usefulness of the 
individual pure states which comprise the superposition. The following definition serves to 
characterize as “metrologically useful” the multimode, equal weight superposition states in 
separable Hilbert space that are extensively more useful for displacement estimation of a 
pre-defined set of self-adjoint generators than the component branches. 
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Definition 1: An equal weight quantum superposition of q linearly independent pure 
states, leu) oc \'^j) ^ is considered metrologically useful when the following 

eondition on the quantity N^^{\u))) is satisfied: 




i Ei=i T^axHeAuioc. (^i I 


e C»(ntot). 


(14) 


where (•|(Ai/)^|-) := = (<^| ® IlAr-i|tu) is the expeeted total 

photon number, A is an algebra of observables on P{C), and Ai-ioc. is the linear subspaee of 
A^^ in which each element is “1-local, ” i.e., has the form E^i ® Ia^-i for Xj G A. 


The set ^i-ioc. should be such that the denominator in Eq. ffT^ is nonzero. The restriction 
to 1-local observables in Definition 1 allows one to use product states as a scaling standard. 
Specifically, given H = ^ -^i-ioc. (here, we have omitted the identity operators for 

clarity), and a product state Ifi) = ® ® it follows that: 


{fiKAHYlfi) < (15) 


Hence, the variance of a measurement of a 1-local observable always scales linearly in the 
number of modes when the system is in a product state. A pure state |a;) having the above 
form is metrologically useful if there exists H (having the 1-local form above) such that 


{uj\{AHfi\oj) e 0{Nntot^BXj^k{'fij ))■ 

The quantity was orig inally introduced as a measure of macroscopicity for quantum 
superpositions in (C^)®'^ |22|; in that context, A = 5u(2,C) (represented by the Pauli ma¬ 
trices) and Utot. is taken to be equal to the number of modes, N. The notion of metrological 
usefulness in Definition 1 refers to the greater ultimate precision achievable in the quan¬ 
tum Cramer-Rao bound when the displacement parameter is encoded in the equal weight 
quantum superposition state |a;) compared to the ultimate precision achievable when the 
displacement parameter is encoded in branches {\'f’m)}m=i comprising |a;). It should be 
noted that one can speak of a superposition state as being metrologically useful only if the 
algebra A is specified. In addition, there may be many ways to write |a;) as an equal weight 
superposition of pure states; in this case. Definition 1 clearly refers to the metrological use¬ 
fulness of |a;) relative to a given decomposition of |a;) into branches. In realistic parameter 
estimation protocols, the branch decomposition could be imposed by the preferred basis of 
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an experiment. 

In the case of a separable Hilbert space the algebra A does not have to be represented 
in the von Neumann algebra B{%). Many observables of interest, e.g., the quadrature 
operators and the photon number operator, are unbounded on T-L but appear in quantum 
optical Hamiltonians of interest to quantum metrology. However, for most quantum optical 
states of interest, these unbounded operators have finite second moment [251 and so have 
bounded variance in these states. In particular, if an unbounded, essentially self-adjoint 
operator x = satisfies {u)\x‘^\u)) < oo for the normalized superposition state |a;) 
then {'il)j\x‘^\ipj) < oo for all j. This feature can be used to introduce a Lie algebra for which 


the state uj is metrologically useful pj(]. A is then formed by taking the 1 -local sums of 
essentially self-adjoint elements of this Lie algebra (we have here assumed a representation 
on %). In order for the denominator of the expression for to be well-defined, at least 
one branch \'ipj) of jcn) must not be an eigenvector of all essentially self-adjoint elements of 
the Lie algebra. 

A simple nontrivial example shows that |ECSj!^(Q;)) is metrologically useful when Ai-ioc. 
is formed from observables of the Lie algebra [)3 = (span{a’^, o, I}, [•, •]), represented as linear 
operators on £^(C) in the usual way. Given a G C, the even and odd coherent states \'ip±) 
(which coincide with |ECS^^;^(«))) exhibit a order |«p variance for measurements of the 
quadrature corresponding to the direction Arg(«) and exhibit squeezing in the variance of 
measurements of the conjugate quadrature corresponding to Arg(a)-|-| |28|. Physically, this 
is due to the fact that the quantized electric field is tt phase-shifted (in expectation) between 
I a) and |—a) coherent states. Explicitly, taking a G M and the quadratures as above, 
then ('i/i+|(Aa;^°^)^|' 0 +) = a^{l + tanha^) -|- 1/2 while = 1/2 — «^(1 — 


tanha^). Thus, if a is purely real, the 6 = 0 quadrature exhibits large fluctuations, while 
the conjugate quadrature fluctuates just below the vacuum level. Since the only observables 
arising from f )3 are the oscillator quadratures and the identity, it is clear that an observable 
E^i exists that exhibits variance on the order of A^a^maxj \zj\‘^ in |ECSjv(a)). 

On the other hand, since every quadrature x^^^ exhibits variance of 1/2 in the coherent 
state |±q;), any 1 -local observable has variance of order A’maxj \zj\^ in 

|±«)®^. Using Definition 1, we see that taking Ai-ioc. fo be composed of observables from 
1)3 allows |ECS^(q;)) to be considered metrologically useful. In particular, by taking Zj = 
X G M, |ECS^(«)) are metrologically useful for estimation of global amplitude displacements 
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0f=i D, (x) of A^-mode oscillators. The estimation of arbitrary local displacements in the 


complex plane comprises a multiparameter {2N real parameters) estimation task j29|]. It 
is an interesting problem whether a measure analogous to can be used to identify 
Schrodinger cat states as a resource for parameter estimation of more general quantum 
dynamics. 

It should be noted that |ECS^(q;)) is not metrologically useful when .Ai.ioc. is composed of 
observables from the oscillator Lie algebra = (span{oio, a\ a, ][},[•,•]) instead of This 
is because the 1-local photon number operator a}jaj exhibits extensive variance A^| 


\a\ 


in the coherent states so that the ratio in Definition 1 exhibits linearly scaling with 

N^ the number of modes, and not the total number of photons. 

We now detail the argument that |HCS^(a;)) are metrologically useful when the algebra 
A is the Lie algebra of observables of 51(2, C). Algebraically, it is simpler to show this fact 
for a closely related hierarchical cat state. By returning to ([2]) and taking |0) = 1?/^+) and 
U = + \'4’-){'4’+\)j the following state is produced: 


|D(a)) = 


1 I ^ |q;) -|- I — Oi) 
\/2 V W2 + 2e-2|“P 


®N 


+ 


\ia) 


-la) 


\/2 - 2e-2|“l' 




(16) 


Consider the 1-local Hamiltonian ® Ijv-i as would describe two- 

photon parametric downconversion into N modes, each with classical pumping amplitude 
z £ C. Because a^\'ip+) = a‘^\'ip+) and clear that 

for any states |^i), |^ 2 ) G spang.{|'i/i+), the following Pauli matrix/two-photon 

quadrature duality holds: 


= + “")!&> (17) 

where ^g appropriate Pauli matrix in 

spang 7 {|'i/i+), From this, it is clear that the variance of ® 

Ijv-i in |D(«)) should be of order A’^| 2 :p|a;|"^. In fact, for a G M the variance is: 

4N‘^'Re{za^)‘^ -(4Re(2^a^) — 8Re(2«^)^ 

Li 

+ 4|2p|aqtaiih c? + cotha^} + 2 | 2 |^) 


( 18 ) 
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which is on the order of A/’^| 2 ;p|a;|^ for Arg^; = 2Arg«. In addition, the variance of 


0 Iat-i in the product states or (e*™ jg most 

of order A^|zp|Q;p, as can easily be verihed. The hnal step in hnding a minimal algebra A 
which allows |fl(a;)) to be metrologically useful is to append the element (l/2)a’^a + 1/4 to 
the set {a^^/2,a^/2} and check that the 1-local observable given by, e.g., ‘ioes 

not exhibit fluctuations in either of the branch states or scaling as 

|«|^. If the 1-local photon number operator were to exhibit such fluctuations, then the ratio 
on the left hand side of Eq. fTlUl would lose the property of scaling with the total expected 
number of photons in |fl(Q;)). It is simple to verify that the 1-local photon number operator 
exhibits variance of order A”lap in these product states and hence, |f2(a)) is metrologically 
useful when A = 51(2, C) := (span£;{(l/2)a’^a -|- 1/4, l2, a^/2}, [•, •]). 

In particular, when 5 G M the calculation above shows that |HCS^(a)) are metrologically 
useful for displacement estimation, where the displacement parameter now corresponds to 
the global squeezing amplitude ^ G M. It is intriguing that while some superpositions of 
5p(A, C) Barut-Girardello coherent states do exhibit squeezing, the |HCS|j(^(a)) state does 
not; in addition, the product states comprising each of the branches of |HCS/^(a)) exhibit 
negligible squeezing if |ap > 1. 

However, if squeezed states and their superpositions are available, one may wonder if there 
exist other types of hierarchical cat states having scaling exponentially in a squeezing 
parameter when the observables of ^i-ioc. are taken from 51(2, C). Indeed, it is known that 
squeezed states provide a higherprecision in the estimation of a single mode squeezing 
parameter than coherent states [SO]. The following hierarchical cat state, having branches 
composed of superpositions of ideal squeezed states, allows for such scaling: 


(D(a) + T>(-a))A(w) 
\/2 I I \/2-\- 2exp(—2a^e^"’) 


®N 


+ 


{D{ia) - D{-ia))S{w) 
\J2 — 2exp(—2a^e^"’) 



(19) 


where we have assumed a,w E M>o and taken S{'w) := as the unitary squeez¬ 

ing operator. For such a and w, the identity D{a)S{w) = S{w)D{ae^) holds and so 
the above hierarchical cat state can be rewritten as S{w)‘^^\fl(ae'^)). The value of 
S\fl(ae'^)) exhibits the same scaling as the value for |r2(ae"’)), i.e., of order 
Na^e^'^, because S{w) acts on 51(2, C) via the adjoint action. As a final remark, we point 
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out that the coherent states ija) are minimum uncertainty states for the Heisenberg uncer¬ 
tainty relation for observables of 1)3 while the even/odd coherent states \'ip±) are minimum 
uncertainty states for the generalized uncertainty relation for observables of 5u(2,C) |3l| . 
This is not a surprising coincidence, as the dehnition of metrological usefulness (Dehnition 
1 ) requires that the maximal uncertainty in the product states comprising the branches of 
a state having the form ([T]) be extensively smaller than the maximal uncertainty in the 
multimode superposition state. 

Thus far, the discussion of Definition 1 has been mainly mathematical. It is useful to 
mention that from a basic physical perspective, the problem of determining the precision of 
optimal estimation of a global real displacement parameter is equivalent to determing the 
energy-time uncertainty in a given quantum state. Therefore Definition 1 can be reinter¬ 
preted from a physical perspective by stating that jcu) is metrologically useful with respect 
to Hamiltonians H G ^i-ioc. if ds maximal decay rate (i.e., minimal time t for which 
becomes distinguishable from |a;)) is extensively greater than the maximal decay rates of 
the branch states i-e., if |a;) is extensively more sensitive to evolution generated 

by H as compared to the branches. In the particular case of |HCS^(q;)) (|ECS^(q;))), one 
can say qualitatively that its metrological usefulness arises simply because its intermode 
quantum coherence causes the squeezing operation (displacement operation) to change it 
more drastically than the product states (lia)®-^) considered independently. 


B. Entanglement entropy 


Because of the orthogonality of the branches, the mode entanglement structure of |HCS^) 
is the same as that of the GHZ states in Hence, |HCS^) exhibits maximal entangle¬ 

ment entropy in the (C^)®^ subspace spanned by {|ej) (D |ej)}jj=i ^2 introduced in Section HI 
In particular, IHCS 2 (a)) exhibits greater mode entanglement than the subset of entangled 
coherent states that cannot be expressed in the form (El). However, the hierarchical photonic 
superpositions are not maximally entangled states of £^(C)®£^(C). The dissipative dynamics 
of the entanglement entropy of entangled coherent states was studied in Refs.jiS. 32 . 33 ]. 

It is known that a nonclassical product state incident on a beam splitter does not neces¬ 


sarily generate entanglement between the output modes [3^ 


destroys the entanglement of a two-mode squeezed state (35|]. It is easy to see that a beam 


. In fact, a 50:50 beam splitter 
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O 0.8 

0.7 



0 0 


FIG. 3: The entanglement entropy of 17(0)|HCS]!^^2(“))) ^ [0.2,3.0] and 6 € [0.1,vr — 0.1]. 

splitter described by the unitary operation B{ 6 ) = e*®/ 2 (“i“ 2 +a 2 “i) acting on two input pho¬ 
tonic modes maps [HCS^(a)) to the state 


1 


v/2 G - e-4|“l 

1 




2 sinh( 2 |a|^) 




( 20 ) 


For moderately large [a], the exponentially decaying term becomes negligible and one is left 
with an entangled coherent state in the output modes of the beam splitter. Hence, the beam 
splitter does not destroy the entanglement of [HCS^(a)) for any values of the transmission 
and reflection amplitudes. The exact entanglement entropy of 5 ( 0 )]HCS 2 («)), calculcated 
as the von Neumann entropy of the reduced density matrix, is shown in FigjHI for a range 
of real a and 9. Except for low-power (a < 1) [HCS^(a)) states, maximum entanglement 
entropy is maintained throughout the range of transmission amplitudes of the beam splitter. 

The quantification of entanglement in terms of an entropic quantity naturally leads 


to questions about its fluctuations. Entanglement fluctuations can 
root variance of a measurement of the entanglement Hamiltonian 


3 e interpreted as the 


36| . i.e., in terms of 


the reduced density matrix pA of a pure state of 77 a ®'Hb it is given by the expression 
I^Se '■= — ^{paHeY where He '■= — \og 2 pA- We show the entanglement fluc¬ 

tuation of i?( 6 *) IHCS 2 (a)) in EiglSl 

We now show that in the Bell basis {]HCS^(«)), ]HCS^(a;)), [ECS^(«)), 0 
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FIG. 4: The entanglement entropy of the image of |HCS]!^^ 2 (‘^)) under independent amplitude 
damping for F = 0.1, a E [1, 2.5] and t E [0,9]. 

I)|ECS^(q;))}, the hierarchical cat states comprise, in some sense, the most stable entan¬ 
glement resource. We consider each mode coupled independently to a zero temperature 
bath of photons, each bath having absorption rate F, with the non-Hamiltonian part of the 
evolution given by 


P = 


^|o,p(i)at] + [ap,a*|. 


( 21 ) 


i=i 


This is the case of (independent) Lindbladian amplitude damping. For an initial state 
p{t = 0) an entangled coherent state, it follows from the well-known solution of the am¬ 


plitude damping master equation [37| that the t ^ oo asymptotic is unentangled. In 
contrast, |HCS^(a;)) maintain nearly maximal entanglement entropy Se throughout the 
non-Hamiltonian evolution as long as jap > 1, as seen in Fig. IH 

The entanglement entropy for this state was calculated from the analytical expression, 
which we omit. It should be noted that for |q;| < 1, the entanglement entropy decays 
with time, but is still substantial for t > T”^. The persistence of entanglement under the 
amplitude damping map can be simply seen by considering a limiting model. By taking the 
low-power a ^ 0 limit, it is clear that |HCS^(q;)) exhibits an inner product of magnitude 1 
with the state 


1 

72 


(| 0 }® 2 + | 1 )®"), 


( 22 ) 


i.e., a superposition of two-mode vacuum and the product Fock state |1) 0 jl). Under the 
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FIG. 5: The entanglement entropy fluctuations of [/(0)|HCS^^2(o^))) « € [0.2,3.0] and 9 € [0.1, vr — 
O.lj. 

amplitude damping map defined in Eq. limi . the state fl^ evolves to 

i(|0>{0| ® |0>{0| g ^ 

+ (e-""|l)(l| + (1 - e-2''‘)|0)(0|) ® {e-"''‘|l>{l| + (1 - e-2''<)|o)(o|)), (23) 

Taking the partial trace to form pi{t), one computes limt^.oo —tr(pi(t) log 2 Pi{t)) = 1. The 
robustness of the entanglement entropy under amplitude damping exhibited for large |q;| 
is proven by considering the |q;| ^ oo asymptotics. An explicit calculation shows that for 
any finite |a;|, limt^.ooAg = 0, whereas fimf_j.oo lini|a|^.oo <S'e = 1. The increased stability 
of the entanglement of hierarchical cat states to local amplitude damping (relative to the 
entanglement of entangled coherent states) makes these states desirable targets for quantum 
optical state engineering and optical quantum communication. 

IV. GENERATION OF |HCS^(a)) 

In this section, we limit ourselves to proposals for experimental generation of the two¬ 
mode hierarchical superposition state |HCS^(q;)) because the main difficulties are already 
present in this case. For all of the proposals we describe, a generalization to V > 2 requires 
the experimenter to overcome a linear (with N) increase in errors associated with imperfect 
implementation of the required unitary operations, in addition to the usual problem of 
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FIG. 6: Quantum circuit diagram for the transformation \'4>-)i 0 |0)2 ® |V’+)3 0 |0)4 
75 (|V’+)i G |V’+)3 + <8) |V’-)3)- 


decoherence due to photon losses. 

If an experimenter has access to arbitrary unitary operations in the sub-Hilbert space 
spanned by |±q;) over a range of amplitudes a, then IHCS 2 (a)) can be readily generated. 
Specifically, one applies the “7r/2” (or “50:50”) beam splitter ?7i2(7r/2) = to 

the product state oc (jv^Q;)! — |—■\/ 2 q;)i) 0 |O )2 to produce the state oc |a;)i|— q ;)2 — |—q;)i|q;) 2 . 
Applying the phase shift produces the Bell state (l/-\/2)(|t^+)|t^_) -|- Ap¬ 

plying = \'ip+){'ip- \ + \'ip-){'ip+\ on mode 2 produces |HCS^(a;)); alternatively, |HCS^(«)) 
is produced (up to a global phase) conditional on the application of annihilation operator 

for preparing entan- 


a 2 to the above Bell state. Along these fines, the method of Ref. j 
gled coherent states by a coherent photon loss may be modified in a simple way to produce 
the following family of states: 


0 (1^+)®^ ± . (24) 

The method is based on the observation that a coherent photon loss can generate photonic 
HCS from a product of single-mode Schrodinger cat states. For example, the above HCS 
state is equivalent (in projective Hilbert space) to (a —5e“*^)(|'0_)a0|'i/i+)ft). The implemen¬ 
tation of the coherent photon loss via a linear quantum optical circuit is shown in Figjbl In 
detail, we append vacuum modes to the tensor product 0 to form the initial state 
|'0_)i 0 10)2 0 |' 0+)3 0 10 ) 4 . Uij{e) := ig ^ 50:50 beam splitter with e <C 1, i.e., 

the beam splitter is highly transmissive for mode i. Applying 1724(7i'/2)e*'^“2“2^^2(e)1734(e) to 
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the initial state produces: 


e*'^ -h 1 

\ ^ 1 * 

acose/i ® \—a — 

\/^ 

sm 6/2 <8 \a- 

Cb 

< 

1 1 

\ 1 ' 

-h acose/i ® |—a- 

a/9 

sine )2 ® \a- 

Cb 

< 

1 t 


— —Qicose)i ( 8 ) a- 

v 2 

sin 6)2 ® \—' 

, V , e*‘^ + l 


— —acose)i( 8 ) a- 

sme )2 <8 —' 

' ' ^/2 



- 1 

Qi4> + 1 


Sin 6)4 (K) |a cos e/3 


a 


V2 

e# + 1 

f.i 4 > _ 1 


sine)4 ® |—a;cose)3 
sin 6)4 ® |q; cos 6)3 


x/2 


sin 6)4 ® —acos 6 ) 3 . 


(25) 


The coherent photon loss is now implemented by photodetection on mode 2 , modeled by 
application of the annihilation operator a 2 . In the final step, one traces over the modes 
2 and 4. In e ^ 0 limit, the “+” state of expression (j24D is produced; if photodetection 
is carried out on mode 4 instead of mode 2, then the state is produced. For initial 

superpositions \'ip±) having large \a\‘^, e must be concomitantly decreased to maintain high 
fidelity of the output state to |HCS^(a;)). The decrease in e necessarily increases noise in the 
photodetection process. In addition, for large |«|^, it is vital to generate the initial product 
state ® |'^+)3 with high fidelity. The next method that we discuss readily satisfies this 
requirement. 

The experimental generation of single-mode photonic Schrodinger cat states \'tp±) via dis¬ 
persive interaction between the monochromatic electromagnetic field and a superconducting 
two-level system Q or Rydberg atom H provides some clues toward feasible methods for 
preparation of photonic HCS. To extend these protocols to the many-mode case, one must 
effectively entangle the field states of spatially separated resonating cavities. For example, it 
has been proposed to generate entangled coherent states by sequential coupling of a Rydberg 


atom to two microwave cavities 


4l| . In general, proposals for creating entangled field states 


involve coupling the field modes to easily controllable, low-dimensional quantum systems. 

A simple scheme for generating IHCS 2 («)) from a tensor product of even coherent states 
\'ip+) ® |' 0 +) is as follows: 


!</>+) ® 1^5+) ^ ^{|</-+> + IV--)) ® \i’+) ^ |HCS+(a)) 


(26) 


where H := l/y/2{ax + is the Hadamard gate in the subspace /C spanned by orthonor- 
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mal basis of even/odd coherent states and CNOT := ® I + \'ip-){'ip- \ ® is the 

conditional operation on the second held mode. To implement the Hadamard operation, 
it is sufficient to generate the following superposition: 


«\i’+) = 


^ (!</>+) + 
Vl + e-2«" 


+ Vl - 


2 Vl - 


, Vl + e-2«" 

\o') H- -j= 

2\/r 


- Vl - e-2«", 




-a). (27) 


Arbitrary superpositions of photonic coherent states |±a;) can be generated by a dispersive 


coupling between a coherent microwave held |a;) anc 
can be prepared in an arbitrary pure state of C 


la) anc 

. a 1 

of 

39] 


a transmon qubit if the transmon qubit 


39j]. In addition, it has been proposed 


to generate parametrically tuning j4^. It is worth noting that i7|V'+} is an eigenvector 
of the operator |q;)(q;| — |—«)(—q;|, which is proportional to the observable corresponding 
to the measurement which optimally detects |«) or |—a) (in the sense of quantum binary 
distinguishability problem with equal a priori probabilities and Bayes’ cost criterion 0 ) 
with maximal probability of success. The pure states H\'ip±) have been studied for their role 
in optimal detection of coherent states |±a}, i.e., the “binary phase shift key” (44 1. 

The CNOT gate in the scheme (l26l) is more difficult to engineer than the Hadamard gate 
because it requires not only a large intramode coherence time for the even and odd coherent 
states, but also a large intermode coherence between two microwave cavities. However, if 
two transmon qubits can be prepared in a maximally entangled (i.e., GHZ) state of (C^)®^ 
and independently coupled to spatially separated photonic modes of microwave cavities via a 
dispersive interaction, this CNOT gate can be implemented. We now provide the details for 
factoring the unitary operator corresponding to the CNOT gate into easily implementable 
unitary operations on the held/qubit and qubit/qubit subsystems. 

First, note that one can factorize the CNOT gate on /C ® /C into the following product of 
local Hadamard gates and conditional gate: 


CNOT = (I ® /7)(|r/^+)(^+| ® I + ® a,)(I ® H). 


(28) 


We have already described the procedure for applying a Hadamard gate to the field via the 
local coupling of the field mode and transmon qubit; hence, we take the initial state to be 
H\'tp+)i ® H\'tp +)2 ® \g)ai ® \g)ai (where we now explicitly include the field mode labels 1, 2 
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H\y/+)i 

\g)al 

\g)a2 

H\y /+)2 


P 

-f- 


a 


X 




p 



FIG. 7: Quantum circuit diagram for the transformation i7|'0+)i ® i7|'0+)2 ® \g)ai ® Id)ai 

|HCS+(a))i,2. 


and the transmon qubit mode labels Ui, 02 ) and show how to implement the conditional 
gate. Let an orthonormal basis for a transmon qubit Hilbert space be taken as {jg'), |e)}. A 
quantum circuit diagram showing our method for indirectly performing the CNOT gate on 
the initial product state is shown in Fig. [71 In this circuit, the first field/qubit operation is 
a TT rotation of qubit oi conditioned on the parity of field mode 1 and is labeled in Fig. |7| 
by the operation with the P superscript. Explicitly, this unitary is given by 






(29) 


A similar conditional transformation has been achieved experimentally in Ref.(39j]. This 
transformation should be followed by a CNOT gate between the qubit modes Ui and a 2 , 
as shown; we assume that this gate is accessible with high fidelity by precise control of the 
qubit-qubit state. At this point, the full normalized state is given by 


{{\'ip+)l ® \'lp+)2 + |'0+)l ® \'tp-)2) ® \g)ai ® 1^) 


a 2 


+ (| V ’->1 ® l '/'+)2 + ® | e ) o , ® | e ) aj ) . 


(30) 


The next step is a tt rotation of field 2 conditioned on the state of qubit 02 . This operation 
has been implemented in the experiment reported in Ref.(l 


Recall from Section Hill that 


the TT phase rotation operator acts like in the subspace /C. Applying again the CNOT 
gate between qubits Ui and a2, followed by the parity-conditioned tt qubit rotation, gives 
the desired CNOT gate in /C 0 /C. Finally, applying the local Hadamard operator 1® H 
produces the output state [HCS^(a)) 0 \g)ai ® \g)a 2 - 

The above method for generating IHCS 2 («)) is not the most efficient possible. It would 
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be favorable to utilize a single qubit or few-level mode which can be sequentially entangled 
with both fields 


V. MORE EXOTIC HIERARCHICAL SUPERPOSITIONS 


The notion of hierarchical cat states can be extended to deeper levels of hierarchy. The 
principal motivation for an analysis of these states comes from the theory of quantum error 


correction, which makes use of encoc 
unwanted decoherence. In Refs 


ed states to strengthen quantum information against 


46l . l47l | a class of “concatenated” GHZ states of the form 


IC-GHZ^yv) := l/V^dGHZ^)®-^ ± |GHZ^)®-^) were introduced as entangled states which 
are relatively stable to local noise compared to the full GHZ state |GHZ^^). An analog of 
the G-GHZ states of can be constructed from HCS in by forming: 


|c-HCSf,j„) - ^ (|HCS+>®" ± IHCS;;)®"). 


(31) 


This state can retain coherence on the scale of N modes even after global coherence on the 
scale of all MN modes has been lost. The C-HCS states are expected to be useful as encoded 
photonic states for continuous variable quantum error correction schemes. Of course, the 
entangled coherent states can be concatenated in a similar way: 


1 


|C-ECSi,„(a)) := ^ (|ECS+(a))®" ± |ECS;i(a)>®") . 


(32) 


It also follows from the basic theory of quantum binary distinguishability that the 
optimal projection-valued measurement for distinguishing from |—has ele¬ 

ments {|G-ECS+^(a))(G-ECS+^(a)|,|G-ECS^^(a))(G-ECS^^(a)|}. The G-HGS and G- 
ECS states are robust quantum resources in the sense that if coherence is lost among the 
M blocks of N single-mode systems, a statistical mixture of A’-mode entangled states re¬ 
mains. To lose all entanglement, the intermode coherence in "H®-^ must subsequently be 
lost. In a higher-order hierarchical cat state, these “shells” of coherence degrade accord¬ 
ing to the strengths of local and nonlocal interactions. It has been suggested to generate 
|G-GHZj{^^) in spin-1/2 chains by application of the 2-local Mplmer-Sorensen unitary gate 
to the AM-mode GHZ state 1/ £ 7 ]. Efficient preparation of hierar¬ 
chically encoded entangled states of (£^(C))®'^^ represents a great challenge for continuous 
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variable quantum information processing. 
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